TUM-HEP-554/04 
UCSD/PTH 04-13 
hep-ph/0409293 
September 2004 



b 



B Xsin- in the MSSM at NNLO 



Christoph Bobeth"'*, Andrzej J. Buras" and Thorsten Ewerth" 

" Physik Department, Technische Universitdt Miinchen, D-85748 Garching, Germany 
Physics Department, University of California at San Diego, La Jolla, CA 92093, USA 



Abstract 

We present the results of the calculation of QCD corrections to the matching conditions 
for the Wilson coefficients of operators mediating the transition b — > sl^l~ in the 
context of the MSSM. Within a scenario with decoupled heavy gluino the calculated 
contributions together with those present already in the literature allow for the first 
time a complete NNLO analysis oi B ^ Xsl~^l~ . We study the impact of the QCD 
corrections and the reduction of renormalization scale dependencies for the dilepton 
invariant mass distribution and the forward-backward asymmetry in the inclusive decay 
B Xsl^l~ restricting the analysis to the "low-s" region and small values of tan p. 
The NNLO calculation allows to decrease the theoretical uncertainties related to the 
renormalization scale dependence below the size of supersymmetric effects in 5 — > 
XJ^l" depending on their magnitude. While it will be difficult to distinguish the 
MSSM expectations for the branching ratio from the Standard Model ones, this can 
become possible in the dilepton invariant mass distribution depending on the MSSM 
parameters and s. In this respect the position of the zero of the forward-backward 
asymmetry sq is even more promising. 
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1 Introduction 



The recent measurements of the branching ratio of the inclusive decay B Xgl^l^ (/ = e, /i) 
of the Belle Collaboration [1] and the BaBar Collaboration [2] are expected to provide an 
important test of the Standard Model (SM) and possible new physics effects at the elec- 
troweak scale. Furthermore they allow for the extraction of informations complementary to 
those from the radiative inclusive decay mode B — > X^'^ which is nowadays well known both 
experimentally and theoretically in the SM and puts non-trivial constraints on parameters 
of models beyond the SM. 

In the discussion of the decay B Xsl^l^ the major theoretical uncertainties arise from 
the non-perturbative nature of intermediate cc states of the decay chain B X^J/tp — > 
XJ~^l~ and analogous higher resonances. These decay channels interfere with the simple 
flavour changing decay mechanism B Xsl^l~ and the dilepton invariant mass distribution 
can be only roughly estimated when the invariant mass of the lepton pair s = = {pi- +pi+y 
is not significantly away from Mj^^ resulting in uncertainties larger than ±20% [3]. For this 
reason the charmonium decays are vetoed explicitly in the experimental analysis [1,2] by 
cuts on the invariant dilepton mass around the masses of the J/ip and ip' resonances. 

A rather precise determination of the dilepton invariant mass spectrum seems to be pos- 
sible once the values of s are restricted to be below or above these resonances. Then the 
calculation can be performed using perturbative methods whereas non-perturbative correc- 
tions can be addressed within the framework of Heavy Quark Expansion (HQE). However, 
contrary to the semileptonic decay B X^^Ji^i and the radiative decay B Xg'-f this 
method is not applicable in the endpoint region of the spectrum as pointed out in [4]. Here 
other approaches have to be used such as for example Heavy Hadron Chiral Perturbation 
Theory (HH^PT) by summing over the kinematically allowed exclusive channels to reliably 
estimate the magnitude of the endpoint decay spectrum. 

At the moment the low-s region, accessible to / = e and /i, is theoretically best understood. 
The non-perturbative corrections (Aqcd/''^6)" to the dilepton invariant mass distribution are 
calculated up to the order n = 3 [4-8] and turn out to be small compared to the leading 
perturbative contribution - however, still involving poorly known matrix elements of the 
Heavy Quark Effective Theory (HQET) for n = 3. Furthermore the effects related to the 
tails of cc resonances in the low-s region of the decay B Xgl^l" were estimated model- 
independently by employing an expansion in inverse powers of the charm quark mass in [9] 
and the size of these AQ(-;j~,/m^ corrections was found to be similar to the size of AQ(-;j~,/m^ 
corrections. Because of the smallness of the non-perturbative corrections in the low-s region, 
the B XJ^l" decay rate is precisely predictable up to about 10% uncertainty. 
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The calculations of the perturbative contribution [10-12] up to the complete next-to- 
leading order (NLO) in QCD [13,14] in the SM had not reached this precision. In a series of 
recent papers the calculation was extended to the next-to-next-to leading order (NNLO) in 
QCD being almost complete up to the missing two-loop matrix element contributions of the 
four quark operators O^-Oq which are expected to be small^. These calculations comprise 

• corrections to the Wilson coefficients [17], 

• the anomalous dimension matrices (ADM) of the renormalization group evolution of 
the Wilson coefficients [18-21]^ 

• the virtual and real corrections to the matrix elements of the relevant operators [22-24]. 

Within the SM the inclusion of NNLO corrections reduces the branching ratios oi B ^ 
Xse~^e~ and B X^jj,'^ ji~ by typically 12% and 20%, respectively [25]. Furthermore uncer- 
tainties due to the dependence on the renormalization scale of the top quark mass ~ Mw 
become reduced from about ±16% to 3% [17] and the inclusion of the NNLO matrix element 
corrections decrease the low energy scale dependence ~ rriij from ±13% to a value about 
±6.5% [22,23]. Furthermore electroweak corrections were found to be a few percent [21] 
removing the scale ambiguity of Oe when going beyond LO. 

Apart from the branching ratio and the dilepton invariant mass distribution, the differ- 
ential forward-backward asymmetry of leptons represents the third interesting observable 
in the decay B X^l^l" . The leading contribution to the forward-backward asymmetry 
arises in the SM at the NLO and thus the inclusion of the NNLO corrections drastically 
reduces the renormalization scale dependence in predictions of this observable. In particular 
it is very sensitive to new physics effects and further, sq = So/ml, the position at which the 
forward-backward asymmetry vanishes provides an important test of the SM [26] . Within 
the SM the inclusion of NNLO corrections in the evaluation of sq leads to a shift of 10% to 
higher values accompanied by a reduction of the uncertainty due to renormalization scale de- 
pendencies in the prediction from typically 15% to 5% [27-29]. The electroweak corrections 
shift So by ±2% [21]. 

Clearly, in view of the improving experimental situation of the ongoing S-physics ded- 
icated experiments, such as the BaBar and Belle experiments, hopefully the experimental 
uncertainties will decrease. Presently it is desirable from the theoretical side to restrict fu- 
ture experimental analysis of the dilepton invariant mass distribution to regions below and 
above the cc-resonances. 

-"^The analogous corrections to B ^ Xgj are 1% [15, 16]. 

^The self-mixing of the four-quark operators Oi-Og is not published yet [20], however the relevant result 
for B — > Xsl~^l~ can be found in [21]. 
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Besides testing the SM, once the experimental accuracy improves, the inclusive decay 
B — i> XJ~^l~ will also allow to constraint models involving new physics scenarios beyond 
the SM. The reliability of such constraints depend crucially on theoretical uncertainties due 
to higher order corrections in the prediction of observables as demonstrated by the SM 
analysis in the case of the importance of NNLO QCD corrections. In the present work we 
report the results of a calculation of QCD corrections to the matching conditions for the 
Wilson coefficients of operators mediating the transition b sl~^l~ in the context of the 
Minimal Supersymmetric Standard Model (MSSM). We chose a scenario in which the down- 
squark mass matrix decomposes into 2x2 matrices for each generation and furthermore a 
heavy decoupled gluino within the MSSM parameter space ensuring the completeness of the 
calculated QCD corrections. 

The scenario of the MSSM that we study here has been introduced in [30] as "Scenario 
B" (see also [31]). It is a generalization of the "Scenario A" which was used in the context 
of the calculation of the QCD corrections to i? — > Xs ^iuu, Bs^d ■, K — ^ tivv and 

Kl —>■ fi~^fi~ in the MSSM [32]^. In this paper QCD corrections to the relevant penguin 
diagrams, the box diagrams and the neutral Higgs penguin diagrams have been calculated. 
While the results of these calculations were ingredients of a NLO analysis of the decays 
considered there, they contribute to -B ^ Xsl~^l~ first at the NNLO level and consequently 
enter our present analysis. Actually as we concentrate on the region tan/3 < 10, only Z'^ 
penguin diagrams and box diagrams calculated in [32] are relevant here. The large tan/3 
region where also neutral Higgs penguins are relevant, will not be considered here. 

Taking into account the results of [32] , the known results for the Wilson coefficients of the 
magnetic penguins, the NNLO corrections to the matrix elements of the relevant operators 
from [22-24] and their three loop anomalous dimensions calculated recently in [18-21], the 
only missing ingredients of a complete NNLO analysis of 5 — Xsl~^l~ in the MSSM are 
the QCD corrections to the Wilson coefficients of the four-quark operators Oi-Oq and the 
semileptonic operator Og. These missing ingredients are calculated here for the first time. 

The main objectives of our paper are then as follows: 

• the calculation of the matching conditions in question at 0{as) which requires the 
evaluation of a large number of two-loop diagrams, 

• the calculation of the dilepton invariant mass distribution in B ^ Xsl~^l^ and of the 
related forward-backward asymmetry in the MSSM at low s and tan /3, 

• the investigation of the renormalization scale dependence of the observables in question 
■^The analytic results given in [32] are applicable also in "Scenario B" considered here. 



3 



and of the impact of the NNLO corrections on these observables in comparison with 
the NLO results, 

• the comparison of the NNLO resuhs in the MSSM with those obtained in the SM, 

• the comparison of the size of MSSM corrections with the theoretical uncertainties in 
the SM. 

The outline of this paper is as follows. In Section 2 we briefly review the elements of the 
MSSM relevant to the scenario with decoupled gluinos. Section 3 summarizes the low-energy 
effective Lagrangian for the transition b — > sl^l~ and the corresponding Wilson coefficients 
including 0{as) corrections in the context of the MSSM. Section 4 presents the formulae 
for the dilepton invariant mass distribution and the forward-backward asymmetry of the 
leptons in the decay B — > XJ~^l~ including all NNLO corrections. The phenomenological 
implications for both observables will be given in Section 5. We summarize and conclude in 
Section 6. Finally the appendices collect the analytical results of the Wilson coefficients. 



2 The considered Scenario of the MSSM 

Let us start by specifying the scenario of the MSSM in which the analytical calculation will 
be performed.^ First, we take the down-squark mass-squared matrix to be flavour diagonal 
so that there are no neutralino contributions to flavour- changing h s transitions, and 
second, we assume the gluino with mass Mg to be much heavier than all other sparticles. 
This ffist assumption corresponds to "Scenario B" described in detail in [30] which was also 
used in [31]. The second assumption leads us to an "effective MSSM" with decoupled gluino 
at the scale Hg ~ Mg [33]. Neglecting all the l/iWg effects, the only modified couplings 
relevant for the NNLO corrections to B ^ Xsl^l~ come from the "chargino - up-squark - 
down-quark" vertex. 
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*For the notation and conventions of mixing matrices and couplings we will adopt see Section 2 of [32] . 
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It is at this scale jig where these couphngs as well as the up-squark masses rriu and mix- 
ing matrices F^^-^ of the "effective MSSM" are determined in the matching with the full 
MSSM. All of them are understood to be MS renormalized quantities in dimensional regu- 
larization. We refrain here from shifting the up-squark masses and mixing matrices of the 
"effective MSSM" into the on-shell scheme in order to avoid the appearance of large loga- 
rithms "ln(yUg/mu)", as can be seen by inspection of (A.l) and (A. 2). Then the next step 
is to integrate out successively all other particles with masses much smaller than Mg and 
much larger than rrit when going to smaller scales using NLO renormalization group (RG) 
equations between all occurring matching scales. In our analysis, however, we integrate out 
all sparticles other than the gluino in one step with the top quark, taking into account the 
LO RG running between fig and fit for up-squark masses and their mixing matrices F^. Due 
to the quartic QCD-interaction of the scalar squarks the LO RG equations of masses and 
mixing matrices are coupled and found to be 



do '^s 
ji—m^ = — 



6=1 



(2.4) 



7 o 6 6 2 

d ,j a., 8 ^rj mi 



with 



e=l c=l "a 



= T^1]^^,T^\ U^,^6 ^ diag(l, 1, 1, -1, -1, -1). (2.6) 



The down-squark mixing matrix F-^ still retains its 2 x 2 block structure after scaling it 
down from fig to fit using LO RG equations, and thus neutralino contributions are absent in 
b s-|- (light particle) decays in LO electroweak interactions at the scale fit- 
So far neither squark masses nor mixing matrices have been measured, and thus in the 
numerical analysis we would like to vary the fundamental parameters of the MSSM rather 
then the squark masses and mixing matrices of the "effective MSSM" . Since the latter are 
determined from the former, when decoupling the gluino in the MS scheme at the scale fig 
such a RG evolution becomes necessary when calculating the two-loop "matrix- elements" 
at the scale fit when decoupling in a second step the heavy SM particles and the remaining 
(apart from the gluino) sparticles. 

3 The Two-Loop Matching Conditions 

The framework of effective theories applied to electroweak decays is a convenient tool to 
resum QCD corrections to all orders using RG methods [34]. As explained in the previous 
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section, the mass hierarchy of the SM and the considered extension - the "effective MSSM" 
- allows for integrating out the heavy degrees of freedom of masses Mheavy > Mjy . The effect 
of the decoupled degrees of freedom will be contained in the Wilson coefficients of the QCD 
and QED gauge invariant low-energy effective theory with five active quark fiavors. 

The effective low-energy Lagrangian relevant to the inclusive decay B — ^ Xsl^l~ resulting 
from the SM and the considered scenario of the MSSM has the following form 

£efr = >^^qcdxqed(m, s, c, b, e, /i, r) 
4Gf 



V2 



AC 

+ ^ E [^^usVu, + + v:,v,, o„ (3.1^ 



with A = {3 ... 10, EOM vanishing, evanescent} numbering the relevant operators Of and 
the corresponding Wilson coefficients . Here Gp is the Fermi constant and furthermore 
we refrain from using unitarity of the CKM matrix. The first term in (3.1) consists of kinetic 
terms of the light particles - the leptons and the five light quark fiavours - as well as their 
QCD and QED interactions while the remaining terms consist of AB = —AS = 1 gauge- 
invariant local operators^ up to dimension 6 built out of those light fields^. The operators 
of entering the effective Lagrangian can be divided into three classes. 
The physical operators are 

O? = {s^,PLT'^Q){QrPLT''b), = {-s^,PLQ){QrPLb), 

o, = {s^.PLT'^b) ^(g7^T'^g), 06 = {sj,,j,,7,,PLT''b)J2iqr'i'''r'T^q), 

Oj = -^m, {sa^-'PR b) F,,, = -m, (s a^'^P^T'^ b) G;,, 

9s 9s 
2 2 

O, = ^is^,PLb)ihn), 010 = ^(s7^Pl&)(^7^5/), (3.2) 

9s 9s 

where Pl,r = |(1 =F 75) are the left- and right-handed chirality projectors, respectively. 
They consist of the current-current operators {Q = {u, c}), the QCD penguin operators 
C3,...,6 {q = {u, d, s, c, b}), the electro- and chromo- magnetic moment type operators Oj^s and 
finally the semileptonic operators 09,10- It should be noted that the above basis of physical 
^The operators conserve flavours other than B and S. 

^The s-quark mass is neglected here, i.e. it is assumed to be negligibly small when compared to nib- 
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operators results from the SM, however in extensions of the SM other physical operators 
could become relevant, too. In the MSSM scenario chosen here this is not the case for low 
values of tan (3 and the SM operator basis suffices. 

In addition to the physical operators several non-physical operators have to be included 
in the matching procedure of the full and effective theories. The so-called EOM vanishing 
operators that vanish by the QCDxQED equation of motion (EOM) of the effective theory 
up to a total derivative can be found in Section 5 of [17]. They appear in intermediate 
steps of the off-shell calculation of the processes 6 — > 57 and b ^ sg and contribute to the 
final results of Wilson coefficients of physical operators when going beyond leading order 
matching. 

The second group of non-physical operators which have to be considered in the match- 
ing procedure are evanescent operators. Evanescent operators vanish algebraically in four 
dimensions, however in D ^ 4 dimensions they are indispensable and contribute to Wilson 
coefficients of physical operators. We use the same convention for the evanescent operators 
as introduced in the evaluation of the anomalous dimensions relevant to b ^ 57, b sg and 
b sl+l- of [18,19]. 

The specific structure of the operators Oi is determined from the requirement that the 
effective theory reproduces the SM AB = —AS = 1 off-shell amplitudes of 6 — » s+ (light 
particles) at the leading order in electroweak gauge couplings and up to O [(external momenta 
and light masses)^ /M^^g^^y]), but to all orders in strong interactions. The same applies to 
the extensions of the SM. For a detailed description of the two-loop matching of photonic 
AB = 1 penguins {b — > 57) in the SM we refer the interested reader to Section 5 of [17]. The 
matching calculation of the supersymmetric contributions is performed analogously. Here in 
addition helpful details can be found in Section 4 of [32]. 

The Wilson coefficients at the matching scale fit can be perturbatively expanded in as{fit) 
as follows 

C? = + ^Cf W + + . . . , Q = {c,t}. (3.3) 

Contributions to order a" to each Wilson coefficient originate from n-loop diagrams which 
follows from the particular convention of powers of the QCD gauge coupling = y^Anas in 
the normalization of the operators O7,...,io in (3.2). 

The result of the matching computation of the Wilson coefficients of the physical operators 
or can be summarized as follows. At the tree-level the only nonzero Wilson coefficient 
is C2 = — 1. At the one- and two-loop level, the only matching condition in the "charm- 
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Figure 1: Two-loop contributions to the function [Ti]l. The wiggly line denotes the W boson. 
Shaded blobs stand for self-energy insertions with up- and down-squarks in the loop. Possible 
left-right and up-down reflected diagrams are not shown. 



sector" which gets contributions from virtual exchange of sparticles (see Fig. 1) is 



CI 



c(2) 



[Ti]i + m - 



7987 



17 

y 



6 



In 



(3.4) 



with fi]^ being the renormahzation scale in the "charm-sector". In the notation of [17], 
[Ti]j = T{x) with X = mf/M^ is the SM top quark contribution. Due to the chosen 
renormahzation prescription the first diagram given in Fig. 1 is completely "renormalized 
away". Thus C2*'^'* is not affected by virtual sparticle exchange. The last two diagrams 
contribute to [Ti]l for which we obtain^ 



a=l n=u,d 



arcsm ■ 



Xg InXo — 16Xa 



208 



(3.5) 



where Xq^ = m'^^/M^, and the definition of the Clausen function Cl2(x) can be found in 
Appendix A. As far as the remaining matching conditions in the "charm-sector" and the 
function [Ti]] are concerned we refer the reader to [17]. 

The one-loop and two-loop matching conditions in the "top-sector" are 



C: 

a 



t(i) 



t{n) 



a 



t{i) 



t{i) 



0, 



0, 
0, 



a 



t{n) 



a 



t{2) 



t{2) 



10 



t{2) 



16 



15' 



''Here we assumed ruq^ > Mw/2 which is clearly fulfilled. 
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(t{n) 
9 



fi(ri) 
10 



(3.6) 



The various functions [XY in (3.6) indicate their origin when matching the h —>■ s+(hght 
particles) Greens functions of the full and effective theory 

• [A]: on-shell part of IPI b ^ s'~f (see Fig. 2), 

• [B'-'-]: b sl^l^ mediated by box-diagrams, 

• [C"]: b sl^l~ mediated by penguin diagrams, 

• [D]: off-shell part of IPI b 57, contributing to b sl^l' (see Fig. 2), 

• [E\. off-shell part of IPI b — > sg, contributing to 6 — sqq (see Fig. 2), 

• [F\. on-shell part of IPI b ^ sg (see Fig. 2), 

• [G]: IPI two-loop diagrams b sqq (see Fig. 3). 

The index n corresponds to the number of loops in the diagrams which can be classified into 
tree- level (n = 0), NLO (n = 1) and NNLO (n = 2) contributions, see also the comment 
below (3.3). Furthermore each function [X]" receives contributions from different virtual 
particle exchange 



The index i corresponds to 

• i = W: "top quark - W boson" loops (SM), 

• i = H: "top quark - charged Higgs boson" loops, 

• i = X- "chargino - up-squark" loops, 

receiving virtual gluon corrections at NNLO and further 

• 2 = 4: "chargino - up-squark" loops including the quartic squark vertex corrections 
proportional to gg. These diagrams contribute only at NNLO. 





(3.7) 



i={W,H,xA} 
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Figure 2: Two-loop contributions to the functions [Dg]}, [E4]} with i = xA- The circle indicates 
the positions where photons or gluons can be emitted, whereas the square (diamond) indicates 
the positions where only photons (gluons) can be emitted. 





Figure 3: Two-loop contributions to the function [G3]} with z = x, 4. Diagrams with crossed 
gluon lines are not shown. 
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Discarding the contributions {H, x, 4} in the sum of (3.7) one recovers the SM results, 
whereas discarding only {x, 4} one obtains the results for Two Higgs Doublet Models (2HDM) 
of type II provided tan (3 is small. 

Explicit expressions for the various functions [X]" can be found in Appendix A. We stress 
that all parameters appearing there are MS renormalized. To obtain the Wilson coefficients 
in terms of on-shell masses and mixing matrices for squarks, the following steps should be 
performed: 

1. Remove the contributions due to strong quartic squark couplings, i.e. the contributions 
with the index i = 4 in the functions [X]". 

2. Make the following shift of the up-squark mass in the contributions with the index 



Observe that this shift involves only the gluonic corrections, since the contributions 
due to strong quartic squark couplings have already been considered in step 1. 

The above two steps are a direct consequence of the application of the full scheme shift from 
the MS to the on-shell scheme given in (A.l) and (A. 2). 

However, using the Wilson coefficients in terms of on-shell quantities one needs of course 
on-shell input parameters. In our approach (see Section 2) we have MS quantities at the 
scale fit, and shifting them to their on-shell values with the help of (A.l) and (A. 2) only 
reproduced our numerical results in the MS scheme if all squark masses are close in size. 
More properly one should integrate out squarks stepwise if their mass splittings are large, 
and then shift to the on-shell scheme at the appropriate scale for each squark. We chose 
to integrate out all squarks at one scale, and hence we refrain from working in the on-shell 
scheme in our numerical analysis. 

To summarize: 



• The contributions [Aj]j, [Fs]l [Cg]j, [Bio]l with t = W,H,x,A and [Dg]'^, [Bg] 



[EaIw, [G^lly, [Ti]j have been calculated previously with the list of references given 
in Appendix A. 

• The contributions [Dg]}, [Bg]}, [E4]}, [G3]} with i = H,x,4: and [Ti]l have been calcu- 
lated here for the first time with the expressions listed in Appendix A. The contri- 
butions with i = H have been calculated already in [35] and very recently the same 
result has been obtained in [36]. 



z = 




(3.8) 
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4 Differential Decay Distributions 



In this Section we provide the formulae of some differential decay distributions of the decay 
B — s> Xgl^l" . These are the dilepton invariant mass spectrum and the differential forward- 
backward asymmetry with respect to the dilepton invariant mass s of the lepton pair. They 
are given in terms of the Wilson coefficients at the low-energy scale /i;, ~ which are 
obtained by solving the RG equation [17-19] and the matrix elements of the operators 
of the low-energy effective theory. At the scale /Xb usually the rescaled operators Oi = 
as/{4:7i) Oi {i = {7, 8, 9, 10}) are used. The corresponding Wilson coefficients are Ci{^ib) = 
(4 7r)/a, 

The method of the HQE is applicable to the inclusive decay B Xgl^l" predicting 
the leading contribution to be the matrix elements of the quark-level transition h sl^l~ 
whereas non-perturbative corrections of the type (Aqcd/"'^;))" can be taken systematically 
into account. However, this method is not applicable over the whole kinematical range of 
s and in this work we will restrict the analysis to the so-called low-s region [17] below the 
cc-resonances. 

The matrix elements of the four-quark operators Oi {i = 1, ... 6) to the process b — > s/+/~ 
are proportional to the tree- level matrix elements of Oj and Og. It has become customary 
to take them into account by the introduction of the effective Wilson coefficient C^^ and 
Cf. The exact expressions for these effective coefficients relevant for the NLO analysis can 
be found in [13, 14, 17], whereas the NNLO corrections are given in [22] for low values of s. 
This involves an expansion in the ratios rric/mi), A/s/rrib and ^fsj^mc. The calculation valid 
for all dilepton invariant masses s can be found in [24]. 

The dilepton invariant mass spectrum with respect to the normalized dilepton invariant 
mass s = s/ml reads 



ds 

X < (1 + 2s) 



^9 



Air 
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^Gl{mr'y\V*V,f 
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'-'10 



48 7r3 
2^ . 2 a, 



(1 



TT 



CU99(S)] +4 



2 

1 + T 

S 



-eff 



[1 



12 Re ( CfCf* ) [1 



2 a, 



71 



t^79(-s)] 



^pBrems,A ^pBremSjB 
+ 7-. + 



ds 



ds 



2 a. 



TT 



■UJjjiS 



(4.1) 



The functions ujij{s) summarize the virtual and real QCD corrections to the matrix elements 
of the operators Oi (z = {7, 9, 10}) [22,27], whereas the terms dT^'^"^'^^ / ds and dT^''''^'''^ / ds 
result from infrared-finite real corrections [23]. In the numerical analysis we follow [21] con- 
cerning the QCD corrections. However, we will not include the higher order QED corrections 
discussed there, but rather use = 1/133 which yields results close to the once obtained 
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including them, as was found in [21]. 

To obtain the hadronic differential decay rate dT{B Xgl^l") / ds within HQE, (Aqcd/ 
rrih)'^ corrections have to be added to the partonic differential decay rate dT{h sl^l~)/ds 
of (4.1) [4-8]. These corrections were calculated up to the order = 3 [8]. In the numerical 
analysis we will only include the corrections n = 2 as the corrections n = 3 involve poorly 
known hadronic matrix elements. We also include the (Aqcd/'^c)^ corrections of [9]. 

The partially integrated branching ratio Br{B Xgl^l") of the low-s region is 

BriB X m-) = BrjB^X) n'-- dr{B^XJ-l-) d_{fl 

with the boundaries chosen to be g^j^ = IGeV^ and q"^^^ = 6GeV^. A very recent exper- 
imental result of Belle for this quantity can be found in the second paper of [1] Br{B 
Xsl^l~) = (1.49 ± 0.50!*^Q28) ^ 10"^ which is in agreement with the BaBar measurements in 
the second paper of [2] both having comparable errors. 

Commonly the semileptonic decay B — > XJu is used as normalization because the factor 
(m^"'*^)^ - the origin of large uncertainties - cancels in the ratio. An alternative was proposed 
in [37] using the charmless semileptonic decays B Xulv and B XcW in the calculation 
of the inclusive decay B Xg'j reducing the uncertainties due to the charm quark mass 
rric present in r(S — XJu). The application of this method to B Xsl^l~ can be found 
in [21,38] and will be used in the numerical analysis. 

The so-called un-normalized forward-backward asymmetry is defined as 

, Br(B^X) d^T(B ^ Xsl+l-) 

Again the normalization is commonly chosen to be the semileptonic decay B XJiy, how- 
ever also the alternative of the combination of the decays B — > X^li^ and B XJu [21] 
to reduce the uncertainties due to can be applied. The so-called normalized forward- 
backward asymmetry is given by the ratio 

. d^r(B ^ XJ+1-) / dViB ^ x,i+i-) 

The numerator at the parton level of the forward-backward asymmetries introduced in (4.3) 
and (4.4) is 



1 dsdz 4g^3 



X 



-3sRe{CfCl^*) (l + ^/9io(s; 



vr 
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-6 Re{CfC^^*) ( 1 + ^ fnois) ) + A^T^' 



(4.5) 



There z = cos 6 and 6 is the angle between the positively charged lepton and the b 
quark in the dilepton center of mass frame. The functions fij{s) summarize virtual and real 
QCD corrections [27,28]. The real QCD corrections are infrared-finite [29] and 

their contribution does not exceed 1% in the SM. In the following they will be neglected. 
As in the case of the dilepton invariant mass spectrum the non-perturbative contributions 
(AqcD/'^fe)" have to be added to pass from the partonic quantity d'^r{b sl^l~) / ds dz to 
the hadronic quantity d^T{B Xgl^l") / ds dz. They can be found in [4,6,8] whereas the 
(Aqcd/''^c)^ corrections are given in [9]. 

The position of the zero of these asymmetries, sq, is of special interest because it is 
sensitive to new physics. It has the value of sq = 0.162±0.010 at NNLO in the SM. However, 
as a quantity comparable with experiments one should consider = mpo- Therefore an 
additional uncertainty due to the 6-quark mass arises. In [21,24] the value of has been 
calculated at the NNLO in the SM yielding q^ = (3.76 ± 0.33) GeV^ and q^ = (3.90 ± 
0.25) GeV^, respectively, depending on the choice of m^. 

5 Phenomenological Implications 

In what follows we will investigate the phenomenological implications of the MSSM cor- 
rections for the branching ratio, the dilepton invariant mass distribution and the forward- 
backward asymmetry. 

5.1 MSSM Parameters and Constraints 

At the present, neither squark masses nor elements of squark mixing matrices have been 
measured, thus it is more appropriate to scan over the fundamental parameters of the MSSM 
Lagrangian in order to investigate the new physics effects. The special scenario of the MSSM 
under consideration has already been described in Section 2. 

These fundamental parameters determine the masses and mixing matrices of the "effec- 
tive MSSM" sparticle spectrum at the scale fig. We would like to remind the reader that 
the MSSM parameters in [30,32] refer to the so-called super-CKM basis [39] of the scalar 
superpartners of the SM fermion sector. The fundamental parameters of the MSSM relevant 
in our numerical analysis are 

• the charged Higgs mass Mh and tan (3 in the Higgs sector, 

• /i and M2 that parametrize the chargino sector. 
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• the gluino mass Mg ^ jJ-g, 

• the soft supersymmetry breaking scalar masses , of left-handed down- and 
right-handed up-squarks, 

• the soft supersymmetry breaking trilinear couplings Au of up-squarks, 

with , and Au assumed to be real and diagonal matrices. Due to the SU (2) gauge 
invariance M^^ is related to Af^^, namely M^^ = VqkmM'^^Vq^-^. Thus the up-squark 
squared mass matrix cannot be decomposed into three 2x2 block-matrices for an arbitrary 
diagonal M'j^^. 

The decoupling of the gluino requires that the masses of all other sparticles should be 
lighter compared to the gluino mass and consequently effects of order Msparticie/^g can be 
neglected. This provides an upper bound on the sparticle spectrum which is chosen to be 
~ 600 GeV. Further lower bounds have to be fulfilled on sparticle masses by direct searches 
from [40] 

• > 94 GeV for the chargino masses, 

• 'fTT'u > 100 GeV for the 2 lightest up-squarks whereas the remaining squarks are required 
to be heavier than 250 GeV. 

Due to the matching of box-diagrams contributing to b ^ sl'^l^ the Wilson coefficients also 
depend on the masses of sneutrinos. As such contributions are rather small we fix their 
masses to be degenerate, with masses between 100 and 300 GeV. Also the down-squarks 
are approximated by a common mass ~ (300 — 500) GeV, as they only appear in the function 
[Ti]i, which effect is negligibly small. 

We have chosen a scenario within the MSSM with values of tan/3 < 10 to avoid the 
appearance of additional operators which are not present in the SM operator basis (3.2). 

A very important constraint on new physics models is the total inclusive branching ratio 
for B — s> Xs7. It has been shown within scenarios of the MSSM [33, 41] that the NLO 
QCD corrections of one-loop diagrams with virtual sparticles can become important and 
comparable to the present experimental uncertainty of Br{B — > ^^7). Also a correlation 
between the B Xsj and the B Xgl^l^ decays is obvious because both involve the 
Wilson coefficient C7. 

The issue of theoretical uncertainties in Br{B ^sl) is not settled yet. Two main 
points arise here. First the choice of the renormalization scheme of the charm quark mass 
rric in the 2-loop matrix elements of the four-quark operators is still a large theoretical 
uncertainty of 11% [37]. It can only be solved by the calculation of NNLO corrections 
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to Br{B Xg'-i) as anticipated in [42]. The second point is concerned with the model- 
dependences entering the results of Br{B X^'-f) measurements when extrapolating to the 
lower end of the photon energy spectrum in the experimental analysis. In [43] a total inclusive 
branching ratio Br{B X^'j) = (3.34 ±0.38) x 10~^ with a photon energy cut Eq > mb/20 
was quoted. A very recent analysis of the Belle Collaboration [44] uses the full inclusive 
spectrum between 1.8 < < 2.8 GeV, without invoking theoretical models of the photon- 
spectrum. The necessity to introduce the photon energy cut in theoretical calculations in 
order to avoid model-dependent experimental results was also raised very recently in [45]. 
The method proposed there results in larger uncertainties of the theoretical prediction of 
the order of 25%. In our numerical analysis the most recent SM calculations [15,37] will be 
used, however with Eq > mfe/20, and the rather conservative interval 

2.0 X 10^^ < Br{B X,7) < 5.0 x 10~^ (5.1) 

to show the correlations with the B XJ'^l' observables. 

The values of the SM parameters are taken to be as in [21] throughout the numerical 
analysis. 

5.2 Results 

We find that the branching ratio Br{B Xsl~^l~) receives only small corrections within the 
considered MSSM scenario. This is illustrated in Fig. 4 where for randomly chosen points of 
the MSSM parameter space, fulfilling the lower sparticle mass bounds, the resulting Br{B 
Xs'~f) versus Br{B Xsl^l~) is shown. The vertical lines correspond to the SM prediction 
of Br{B Xgl^l^) and the corresponding estimate of the theoretical uncertainty [21]. The 
horizontal lines indicate the SM prediction and theoretical uncertainties of Br{B — > X^'^) 
[15,37]. Deviations are possible from the SM central value Br{B Xsl~^l~) ~ 1.60 x 10~^ 
up to ±(15 — 20)% respecting the experimental bound from Br{B Xg'-y). Therefore the 
observable Br{B XJ~^l~) of the low-s region will not serve as a good candidate allowing 
to distinguish the SM and the considered MSSM scenario in view of the present theoretical 
uncertainties. The reason is the smallness of the MSSM contributions to Cg^ and C^q which 
dominate in the expression for Br{B XJ~^l~) in the low-s region. Although Cf could 
receive a larger MSSM contribution its magnitude is strongly constraint by the measured 
value of Br{B Xg'-f).^ Furthermore, the contribution to [C^'^^l^ to the differential branching 
ratio falls like 1/s and therefore only dominates for values of s < 0.05 which coincides with 

*It should be stressed that this is a quite loose terminology since for the LO expression of Br{B Xg'y) 
the initial Wilson coefficients of the two operators O7 and Os enter. At the NLO this becomes even more 
involved. For a model-independent analysis of this subject in the presence of new (scalar) operators see [46]. 
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Figure 4: Br{B Xg'j) versus Br{B Xgl^l") for randomly chosen points in the parameter 
space of the MSSM scenario. The three vertical lines indicate the SM prediction of Br{B 
Xsl'^l~) [21] and the three horizontal lines the one for Br(B Xs'j) [15,37]. 



the lower end of our integration range. The interplay between various contribution to the 
differential branching ratio within the SM is depicted in Fig. 7. There also a specific point 
in the space of supersymmetric parameters with significant corrections to C^q is shown. 

The position of the zero of the forward-backward asymmetry sq represents a more sensitive 
observable than Br{B — Xgl'^l^) in the considered MSSM scenario. In Fig. 5 we plot 
\/Br[B — s> Xs7) versus Sq of the normalized Ay-b{s) for randomly chosen points of the 
MSSM parameter space. There the vertical lines correspond to the SM prediction of Sq 
and its uncertainties [21,24] and the horizontal lines as in Fig. 4 to the SM prediction of 
Br{B ^ Xsi). 

We note that the points in both plots in Fig. 5 are clustered along a straight line, exhibiting 
very clearly the correlation between the value of Br{B Xg'y) and Sq within models with 
minimal flavour violation (MFV) as pointed out in [47]. 

The straight lines in Fig. 5 are to a very good approximation model independent within 
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Figure 5: ^/BriB~^'J(/y) versus sq, the position of zero of the normalized Afb{s) for randomly 
chosen points in the parameter space of the MSSM scenario. The three vertical lines indicate the 
SM prediction of Sq [21,24] and the three horizontal lines the one for Br(B Xg'y) [15,37]. 
In the left plot the MS charm quark mass is used for Br{B — > ^^7), whereas in the right plot 
mP°'°, resulting in a smaller prediction. 

the class of models with MFV. Only different points on them correspond to different models 
and/or different sets of parameters in a given model. On the other hand the position of these 
lines depends on the parameters of the low energy theory, in particular on the charm quark 
mass that enters sensitively the evaluation of Br{B Xg'y) [37] but is practically irrelevant 
for §0- In the left plot in Fig. 5 we used the MS mass mc(mc) and in the right plot the 
jj^poie j-Qgi^gg^ that results in a different straight line. The SM prediction for Br{B ^^7) is 
lower in the right plot than in the left plot. It is clear that the usefulness of the correlation 
between the values of Br{B ^s7) and sq in testing the MSSM will depend on the progress 
in NNLO calculations for B Xg'j that should significantly decrease the sensitivity due to 
the choice of rric. 

As seen in Fig. 5, in addition to dense points in the ballpark of SM expectations, there 
are values of Br{B — > Xg'j) and sq within the MSSM that are larger and smaller than the 
SM predictions. This should be contrasted with the result in a model with one universal 
extra dimension in which only smaller values of Br{B ^^7) and sq were possible [47]. 

In Fig. 6 we show Br{B Xgl^l^) versus sq. In the left plot the MS definition was 
used for the charm quark mass in the evaluation of Br{B ^s7) whereas in the right plot 
the pole-mass definition. As a consequence the allowed range of the position of the zero of 
Ay-b{s) becomes shifted a bit towards higher values. The comparison of Fig. 5 and 6 shows 
that the position of sq is much more sensitive to the Wilson coefficient C7 and consequently 
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Figure 6: Br{B XJ~^l~) versus sq, the position of zero of the normalized ^fb(s) for randomly 
chosen points in the parameter space of the MSSM scenario. The three vertical lines indicate the 
SM prediction of Sq [21,24] and the three horizontal lines the one for Br{B — >• XJ^l') [21]. 



to Br{B Xs-f) than to Br{B XJ+l') itself. 



"PI" 


Mh = 440.11 GeV, tan/5 = 5.01, fi = -122.87 GeV, M2 = 184.56 GeV, 
Au = diag(370.29, 79.60, 535.71) GeV, 
Md^ = diag(299.63, 127.20, 454.43) GeV, 
Muj, = diag(219.96, 519.91, 167.68) GeV 


"P2" 


Mh = 248.34 GeV, tan/3 = 2.56, /i = 192.83 GeV, M2 = 489.68 GeV, 
Au = diag(-419.30, 525.64, -540.81) GeV, 
M^^ = diag(339.09, 128.18, 393.52) GeV, 
Muj, = diag(232.08, 351.41, 234.77) GeV 


"P3" 


Mh = 451.74 GeV, tan/5 = 4.89, /i = -540.06 GeV, M2 = 582.50 GeV, 
Au = diag(-375.95, -324.59, -497.23) GeV, 
Md^ = diag(503.97, 281.42, 264.06) GeV, 
Mu^ = diag(444.06, 186.86, 417.40) GeV. 



Table 1: Three selected points. In all points Mg = 1 TeV, and the down-squarks and sneutrinos 
are assumed to have masses about = 300 GeV and = 250 GeV. 

In Fig. 7 we show the four main contributions due to \C^%io\^ and Re (cfCf*^ to the 
differential branching ratio, see (4.1), as functions of s for the fixed MSSM parameter point 
"PI" defined in Table 1. Each plot shows the SM (light grey band) and the MSSM contribu- 
tion. To demonstrate the reduction of the renormalization scale dependence fit we show the 
MSSM result when including all calculated corrections (dark grey band - "MSSM") and the 
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Figure 7: Various contributions to the differential branching ratio as functions of s. 



partial MSSM result (shaded bend - "mssm") obtained by discarding all contributions with 
n = 2 and i = {H,x,4} to the functions [X]" in (3.7), but not to the SM. The bands are 
obtained by varying the renormalization scale fit G [120, 300] GeV and the low-energy scale 
fib £ [2.5, 10] GeV. Large deviations from the SM appear in the contribution jCJ'o ^ mainly 
due to the Z°-penguin function [Cg] which is suppressed in |C|^p as can be seen in (3.6). 
The inclusion of the NNLO matching conditions in the MSSM reduces the renormalization 
scale dependence to comparable size as obtained in the SM calculation. 

The sum of this four separate contributions (and the bremsstrahlung contributions) adds 
up to the final differential branching ratio shown in Fig. 8 in the left plot. As before the 
bands are obtained by variation of the renormalization scales fit and fib- The reduction due 
to MSSM contributions is roughly —30% for values of s > 0.15 as can be seen in the right 
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Figure 8: The differential branching ratio for fixed MSSM parameter point "PI" compared to the 
SM result and the partial MSSM result as a function of s (left plot). The relative size compared 
to the SM is given in the right plot. 



plot of Fig. 8 where the relative size compared to the SM result (obtained for fit = 120 GeV 
and fib = 5 GeV) is given by the quantity AQ = Q/Qsm — 1- Thus the shape and magnitude 
of the dilepton invariant mass distribution provides in certain regions of s a more sensitive 
observable then the integrated branching ratio itself in the search for deviations from the 
SM prediction, depending on the MSSM parameter point. It should be noted that the very 
small scale dependence around values of s ~ 0.05 are due to accidental cancellations between 
the 4 separate contributions in (4.1). 

In Fig. 9 we show the normalized forward-backward asymmetry ^fb('5) for the low-s 
region. The left plot illustrates the result for the fixed MSSM-parameter point "P2" and 
the right plot for "P3" that are given in Table 1. The SM result is shown in both plots 
for comparison. Again the bands are obtained by varying the renormalization scales fit and 
fib as in Figs. 7 and 8. Due to the strong correlation of the position of the zero sq and 
Br{B Xs'~f) in the considered MSSM-scenarios further shifts to the left or right (as shown 
in the two plots) of Sq are unlikely. 

In Fig. 10 the fundamental MSSM parameters fi and [Af/jaa are shown versus the position 
of the zero of Afb{s), Sq, for the sample of random MSSM points given in Fig. 5. The lower 
and upper bounds of Sq present in both plots are evidently due to the strong correlation to 
Br{B Xs'^). The "hole" in the fi distribution for values < 100 GeV comes of course 
from the bound on the lightest chargino mass. As can be seen for small values of Sq also 
smaller values of fi are preferred. The allowed values of [Aj/Jsa versus sq generated during 
our random scan are shown in the right plot. Almost no bounds are found here, only towards 
smaller values of sq very small values of [A^/jas seem to be excluded. We could not find such 
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Figure 9: Normalized forward-backward asymmetry AfB{s) versus s in the low-s region for two 
(left, right) fixed MSSM-parameter points "P2" and "P3" (see Table 1) compared with the SM 
prediction. 



correlations for all other soft-SUSY breaking parameters. 
6 Summary and Conclusions 

In this paper we have presented for the first time complete NNLO QCD corrections to 
B — i> Xgl^l^ in the context of the MSSM within a scenario as defined in Section 2. We 
have calculated the missing ingredients of a complete NNLO result and including also con- 
tributions present already in the literature, we were able to calculate with this accuracy 
the branching ratio for B Xgl^l" in the low-s region, the corresponding dilepton invari- 
ant mass distribution and the forward-backward asymmetry. The presented results can be 
applied to all MSSM scenarios with a flavour- diagonal down-squark mass-squared matrix, 
as long as the gluino is heavier compared to the remaining sparticle spectrum and tan/3 is 
small. 

This calculation was motivated by the fact that in the SM the B — > Xgl^l'^ observables 
suffer from sizable renormalization scale uncertainties which are reduced considerably at 
NNLO. Consequently in order to have a chance to see supersymmetric effects of in this 
decay, it is essential to reduce renormalization scale uncertainties in the MSSM as well. 

The results for the relevant Wilson coefficients are collected in Section 3 and the Appendix 
A, where we stated explicitly which corrections were already calculated previously and which 
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Figure 10: Correlation between sq and parameters fi and [A^/Jas produced in our scan over MSSM 
input parameters. The distributions for the other soft-breaking parameters are flat. 



are new. The numerical analysis of the quantities of interest is performed in Section 5. Our 
main finding's follows: 

• The /i dependence present in all quantities of interest at NLO is visibly reduced at 
NNLO depending on the magnitude of the MSSM contribution for the particular MSSM 
parameter point and it is typically of the same size as the one of the corresponding 
SM result. 

• Supersymmetric effects in the branching ratio amount only to at most 20% and conse- 
quently in view of theoretical uncertainties in this quantity it will be very difficult to see 
them unless experimental and theoretical uncertainties will be significantly reduced. 
In this respect the dilepton invariant mass distribution can offer in certain regions of s 
the possibility to distinguish the supersymmetric effects from the SM prediction. Such 
effects can reach up to 30% depending on the MSSM parameters and the value of s. 

• The best chance to observe supersymmetric effects in this decays is through the 
forward-backward asymmetry. We find that the position of the zero Sq in this asym- 
metry can be significantly shifted both downwards and upwards relatively to the SM 
expectation. These shifts are accompanied by the shifts in Br{B — > Xg'-f) as pointed 
out in [47] and shown in Fig. 5. As the predictions for sq is theoretically rather clean, 
accurate measurements could be able to detect possible departure of its value from the 
SM prediction one day. 



23 



Acknowledgments 

We thank Jorg Urban, Michael Spranger, Janusz Rosiek and Mikolaj Misiak for discussions. 
CB. and T.E. have been supported by the German-Israeh Foundation under the contract 
G-698-22. 7/2002. Further C.B. has been supported by the Department of Energy under 
Grant DE-FG03-97ER40546. The work presented here was also supported in part by the 
German Bundesministerium fiir Bildung und Forschung under the contract 05HT4WOA/3 
and DFG Project Bu. 706/1-2. 

Appendix A Wilson Coefficients 

This appendix summarizes the matching results relevant for B Xsl^l^ in the SM and the 
considered scenario of the MSSM as introduced in Section 2. It provides the formulae for 
the functions [X]" introduced in (3.7). 

Dimensional regularization with fully anticommuting 75 and the MS scheme was used 
for all QCD counterterms, both in the full and in the effective theory for light degrees of 
freedom. The only exceptions were the top quark and squark loop contributions to the 
renormalization of the light-quark and gluon wave functions on the full theory side. The 
corresponding terms in the propagators were subtracted in the MOM scheme at = 0. In 
consequence, no top quark and squark loop contribution remained in the "light quark - W 
boson" effective vertex after renormalization. 

The only relevant off-shell electroweak counterterm in the full theory proportional to spb 
was taken in the MOM scheme as well, at = for the s<^b, and at vanishing external 
momenta for terms containing gauge bosons. 

As a consequence of this special choice of the renormalization, all masses of quarks and 
squarks, as well as the mixing matrix and the effective couplings Xl^^ and X^^^ appearing 
in this Appendix are MS quantities. The masses of particles which do not interact strongly 
are not renormalized and thus might be interpreted as their tree-level masses. 

As already stated in Section 3 all the functions [X]\ are equal to zero choosing an on-shell 
renormalization prescription for squark fields and masses [48] and the mixing matrix F^ [49]. 
For example this can be seen by means of the following transformation formulae between 
MS and on-shell scheme, 





(A.l) 
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On the left hand side of these equations the squark masses and the mixing matrix are 
running MS parameters, whereas on the right hand side they take their on-shell values. We 
note that the couplings Qg and ay given in (2.3) are still MS renormalized working with 
on-shell squark parameters. 

We define the mass ratios 

ml _ ml _ Ml _ m\^ _ mj^ 

^~m|;' ^~m|:' ""''"Mf' y^^'Mf/ ""^^'W/ 

and introduce the abbreviations 

'"'-'4: -w^: 

In these equations mt denotes the top quark mass, m^^ and m^^ up and down squark masses, 
Mvi/ the W boson mass, the charged Higgs mass, M^. the chargino masses and finally 
m,>^. sneutrino masses. 

The integral representations for the functions \ji2{z) and Cl2(a;) are as follows 

U,{z) = - r dt^^^^, (A.5) 

Cl2(x) = Im [Li2(e")] = - / rf^ln |2sin(^/2)|. (A.6) 

Jo 

The calculation was performed in the background field formalism in an arbitrary i?^-gauge 
for the gluon gauge parameter and in the t'Hooft-Feynman gauge for the W boson gauge 
parameter. 

Appendix A.l i = W — "top quark — W boson" 

The evaluation of Feynman diagrams contributing to 6 ^ s+ (light particles) Greens func- 
tions within the SM mediated by "top quark - W boson" loops yields the functions denoted 
by the index i = W in (3.7). The explicit form can be found in [17] by using the equalities 

(^), [F,]w = ^o(^), 

B{{x,-\l [Clt^ = Ci{x\ (A.7) 
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The functions [X]]^ have been first calculated in the following papers: [Ay]]^ and [-Fg]^ 
in [33,50-53], [B'^Yw and [C^Yw in [54-57] and [D,Yw, [E^]^ and [^3]^^ in [17]. 



Appendix A. 2 i = H — "top quark — charged Higgs" 

The evaluation of Feynman diagrams contributing to 6 ^ s+ (light particles) Greens func- 
tions within the MSSM (but also 2HDM of type II) mediated by "top quark - charged Higgs 
boson" loops and denoted by the index z = if in (3.7) yields 
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6(s/-l)4 



4(3/-l)3 



7y^^l6y-'+81^/ 
3{y-l)3 



I „„^2 q f -13y4+l7^3_3o^2 / l\ -468y4+321y3_2i55^2_2^ 
+ COt 3{^i ^^^y~y)^ 108{y-l)5 



4451^"^ -7650</^+18153^/^ - 1130^ 
648(s/-l)4 



17y3-31;/2 1 , 7y4-18y3 +261^2+38;; 
SOPTp— y 18(y-l)4 



(A.19) 



30y^~my^-^&y 
81(y-l)4 



Gy^'- 187;/^ +213;/ 
81(y-l)3 



+ 



-30y2+20y 
27(y-l)4 



-35y''+145y^-80y 



(A.20) 



The following functions [X]\j have been calculated previously: [A-j\\j and [-Fs]// in [33,52,58] 
and [-Bjol/f ^^"^ [^gl/i' i'^ [32]- The function [-Dg]^ has been calculated in [35] and confirmed 
in [36]. The results for the functions [E^\j and [6*3]}^ are new. Note that [-B9])/ and [-B'oIh 
vanish due to the approximation of vanishing lepton masses. 



Appendix A. 3 i = X ~ "chargino — up squark" 

The evaluation of Feynman diagrams contributing to 6 ^ s+ (light particles) Greens func- 
tions within the MSSM mediated by "chargino - up squark" loops and denoted by the index 
i = X in (3.7) yields 



[^4 = -EE # hf\y.^ + ^[X^^]2.[X^].3 hf\yA 

i=l a=l Xi ^ J 



(A.21) 



rr^ll nO ^ \^ ST t^f ^^laJ^*^^ J'^S 

= ^''^ 



i,j=l a=l 6=1 Xi 
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(A.22) 
(A.23) 



2 6 , ,2 



m = -l^l^ j^l^i ha[XnaS h^iVa 
1=1 a=l 



(A.24) 



i=l a=l Xi 



(A.25) 



i=l a=l Xi ^ 



5 \yai] 



%.r^C/itn r-^URl ,(0)/ 



nib 



(A.26) 



2 6 



MS, 



i=l a=l Xi 

X ([Xf-^]2.[Xf-].3 /^S'Hya.,L,J + ^[xf-t]2.[Xf«]„3 h^\ya,, L^J } , (A.27) 
I rrib 



[^9,10 



i,j=l a=l 6=1 



f8^\xji, yai, Vbi) +4(1+ yai^— ] fl^'^Xji, yai, Vbi) 



d 



dya 



fg^\xji, Vai, Vbi) +4(1 + VaiW— ] f&\Xji, yai, ^bi) 



d 



dya 



i-(O)/ 



(A.28) 
(A.29) 



• 1 ^K 

1=1 a=l Xi 



(A.30) 



1=1 a=l Xi 



(A.31) 



K-EE^ 

i=l a=l Xt 
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X ([Xf-']2a[Xf^]a3 h^\y^,, + ^[X^^^'uxf-U h^\ya^,L,j] , (A.32) 



YWj^ f;^ti r^C/Li uWl 



= '^EE M^t^^"" ]24^fla3 4'Hya.,L,J. (A.33) 
i=l a=l Xi 

The following functions [X]^ have been calculated previously: [^7]^ and [Fs]^ in [33,41] 
and [B^^oW and [Cf]^ in [32]. The resuhs for the functions [B^J]\, [Dg]^, [^4]^ and [Gg]^ are 
neti;. The expressions for the functions [C)^'']^ and [C^'']^ correspond to the leading and the 
next-to leading contributions to the function [Cf]^ = [Q^]° +«s/(4vr)[Q^]l given in (3.14) 
of [32]. 



Appendix A. 4 i = 4 — "chargino — up squark (quartic)" 

The evaluation of Feynman diagrams contributing to 6 ^ s+ (light particles) Greens func- 
tions within the MSSM mediated by "chargino - up squark" loops containing the quartic 
squark vertex^ instead of gluon corrections and denoted by the index z = 4 in (3.7) yields 

i=l a,b,c=l Xi 



X |[X^^]2jX^],3 [-q['\ya^,yc^) + ^g^'H^/a., Z/c.)] 



+ MM[xI'l^]^,[XI'%, [-qi'\ya^,yc^) + ^ (l/a. 2/«)] ) , (A.34) 

4 ^ ^ 

[^9, 10] 4 = =^^3 E E E PabybiPbc (1 + ^iij [-^7'' ]2a[Xf'']ci 

« J=l /=1 a,b,c=l Xi 



X 



ya^, ya, W/i)[Xf^^],/[Xf^]/i 
TV^^fio\^^^.ya^,yc^,Vf.)[X^-\f[X;'-]J^ , (A.35) 



i^Si = ^ E E (1 + ^.J [^7"']2d[X 



6 

i,j=l a,...,e,g,k=l 



X <! 2y/xJ'iff\xji,yai,ydi)UjlU*^5ae5gdhl5cl - lf\xji,yai,ydi)V*^Vii6aeSgdSblScl 



^Strictly speaking these matching contributions originate from the part of the quartic squark vertex 
proportional to the strong couphng constant 
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(A.36) 

[^9]; = E ^ P^,yb^P|;^ (1 + L,,) [Xf-^]24Xf^]c3 q^\ya^,yci}, (A.37) 

i=l a,b,c=l 

ml = ^Yl E ^Paby'>^P^cil + L,,)[Xl''%a[Xr^Uql^\ya.,yc^, (A.38) 



j=l a,fe,c=l 

X ([X^^]2jX^].3 g?Hl/a.,yc.) + ^[Xf-^]24X^]e3 ^ (l/a. , Z/c.) | , (A.39) 

ml = 0. (A.40) 

The following functions [X]l have been calculated previously: [^7]^ and [^8)4 in [33] and 
[5fo]^ and [Cf]^ in [32]. The result for the functions [B^^]l, [Dg]l, [E^l and [G^ll are new. 



Appendix B Auxiliary functions 

Here we present explicit formulae for the loop functions hf'\x), hl^\x) and qi^\x,y) intro- 
duced in Appendix A. They read 

<^(-) = lg^ln- + -i^, (B.41) 

^2°H^) = ^f£^lnx+3g^, (B.42) 

hPi^) = -%t%-' lr,x + '%^_'^t'' , (B.43) 

hf^ (^) = In X + -0^, (B.45) 
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(B.46) 



a(I)/'^ „,\ _ -48a:^-104a:^+64x T ; /l 1 ^ _L -378x3-1566x^+850x+86 i„ ™ 

y^^y) - — — ^^H^ ~ x) + 8i(x-i)5 

-56x3+258x2+24x-82 



+ 



9(x-iy 

2060x3 +3798x^ -2664x- 170 
243(x-l)4 

224x^-96x 
9(x-l)3 



81(a;-l)5 
12x3-124x2+64a; ^ 



Li2(l-i) + 



9{x-l)5 
24x3 +352^2 -128x-- 32 



27(x-l)4 



9{x-l)4 



]„ „ I -340x^+132x+40 
^"•^"^ 9(x-l)3 



~24x3+176x^-80x 1„ ™ _i_ -28x^-108x+64 
9(x-l)4 ^^^-^ 9(x-l)3 



y, 



81(x-l) 

- 304x3 _ i36g6x^ +29076X- 12062 
729(x-l)4 



- 108x"t+1058x3 -898x^- 1098X+710 
81(x-l)5 



+ 



540x3 - 972x^ +232X+56 1„ ^ _|_ -664x3+54x^+1944x-902 
81(x-l)5 illX-^ 243(x-l)4 



-562x3+1101x^-420x+101 

54{x-l)'* 



Li2(l-i) + 



- 562x3 +1604x^ -799X+429 
54{x-l)5 



I 17470x3-47217x^+31098x-13447 , 
972(x-l)4 "I" 



89X+55 1„ ^ _|_ 38x3-135x^+54x-821 
27(x-l)5 "^-^ 162{x-l)'' 



4^H^,y) = ng-iy^ Li,(i-i) + 



9x3+46x^+49x 



1\ , 81x3+594x^+1270x+71 



I -923x3 -3042x^ -6921x-1210 , 



54(x-l)^ 

3x 

3{x-l)5 



Inx 



324(x-l)4 

32x2-24x-|^-^j^]^ _ _|_ -52x2-109x-7 



10x^+38x 2; _|_ -7x3+30x^-141x-26 



9{x-l)4 



3{x-l)3 



-20x^-52x 1„ ^ I -2x^+60x+14 
3(x-l)4 ^^^X-f 3(a,_i)3 



3(x-l)4 



i„ ™ _l_ 95x^+180x+61 
"^•^"^ 6(x-l)3 



(X, y) 



-20x3+60x^-60x-20 
27{x-l)4 



-60x^+240x+4 
81(x-l)* 



132x"'-382x+186 
81(x-l)3 



20 In T -I- -20x^+70x-110 
27{x-l)4 ^^^-^ 81(x-l)3 



qi\x,y) 

(l?{x,y) 



qi\x,y) = 



4 


x^ In X 


y2 iny 


3(x- 


y) 


L(X-1)4 




4 




X Inx 


y\ny 


3(x- 


y) 


L(x-i)4 




8 




—x^ Inx 

1 ..^ -1 ^ M 





I 4x^^y^+10xi/^-2y^+10x^y-44xy+10y-2x^+10x+4 
9(x-l)3(y-l)3 ' 



+ 



-2x^y^+10xy^+4i;^+10x^i;-20xy-14i/+4x^-14x+22 
9(x-l)3(j/-l)3 ) 



-12xi/+4y+4x+4 
3(x-l)2(y-l)2 ' 
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x-y) 



—xlnx I ylny 

(X-1)3 ^ (y-l)3 



+ 



-'lxy—4y-4x+12 
3(x-l)2(y-l)2 ' 



:B.57) 



+ 



(6x^-9x2+2) Inx _ (6^3-9^24.2) iny 
(x-l)4 (y-\Y 



104x2y2_2023;,;2_(.gg^2_202a;2^+380xy-154y+86x2-154x+56 



81(x-l)3(.y_l)3 



:B.581 



= 9( 



X-J/) 



In X 



4x2y2 _ i43;y2 -|-22y2 _ 14x2y+52xy-62y+22x2 - 62x+52 



27(x-l)3(y-l)3 



:B.59) 



The functions /g""*, /g""*, /i^o\ /g^^ and /g^' can be found in Appendix B of [32] 



(1) 
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